Theory for Quantum Spin and Vortex Glasses 
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We develop a theory for the quantum vortex glass, with both the coupling strengths and the site 
energies disordered. This model is closely related to XY spin glasses and bosons in random media. 
For properly chosen distributions of the site disorder the phase diagram consists of a superfluid 
phase, and Weak and Strong Glass regions, dominated by long range and local fluctuations, re- 
spectively. The most direct experimental manifestation of the Strong Glass transition is the lack of 
divergence of the corresponding non-linear susceptibility. 
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A particularly active field of research these days is the 
quest to understand quantum phase transitions in the 
presence of disorder at T=0. Theoretical interest fo- 
cussed on one dimensional transverse field Ising problems 
[ 1 1 , on quantum critical behaviour of several spin models 
[2-H], and bosons in disordered media p|. In one dimen- 
sion scaling techniques delivered several exact results. In 
higher dimensions, after the original formulation re- 
cently quantum fluctuations were incorporated as well 
|,§. While it was possible to describe the weakly disor- 
dered phases of the models satisfactorily, our understand- 
ing of strongly disordered regions is still incomplete. Ex- 
perimental studies intensified after the discovery of ma- 
terials which are credible realizations of the spin-glass 
models, such as the LiHoYF compounds ||. In this 
material at the low T glass transition, presumably dom- 
inated by quantum effects, the non-linear susceptibility 
Xni did not diverge, in contradiction to recent numerical 
studies @. 

In the present paper we develop a new technique to 
capture the physics of strongdisorder. In agreement with 
some previous approaches |4| we find two types of glass 
transitions, but this new method enables us to charac- 
terize the critical behaviour even in the case of strong 
disorder. One of our key results is that the non-linear 
susceptibilty does not diverge at criticality. As this re- 
sult is due to local fluctuations, we expect it to hold in 
a wide variety of glasses, Ising as well as the XY type, 
considered here. 

Specifically we consider the quantum gauge glass 
model, which describes hard core bosons propagating in 
a magnetic field, and as such, is believed to contain the 
essential physics of vortices in Type // superconductors. 
H Elastic vortex models in the classical limit are also 
argued to support two distinct glassy phases || , with re- 
cent experiments indicating a transition between them as 
the magnetic field is tuned. 

We consider the Hamiltonian: 



(1) 



furthermore < >= Jq/N, where N is the number of 
lattice sites, and < \Jij\ 2 >= J/N. With this conven- 
tion we not only have random phases along the bonds, 
representing the random fluxes, but also random bond 
strengths. As our results depend only on the second cu- 
mulant of the Ji/s, the two models are expected to give 
the same results, fii = fi + hi, where fi controls the den- 
sity of the bosons, and hi is a random field or site energy, 
with distribution P(h) over the finite support [—A, A]. 
Finally ai and a\ annihilate and create hard core bosons, 
or equivalently correspond to spin-lowering and raising 
operators. 

At sufficiently high values of /z there will be exactly 
one boson per site, leading to the formation of a "Mott- 
insulator" jj| , which is characterized by a gap above the 
ground state and a lack of an order parameter. This 
corresponds to the paramagnet in the spin language and 
to the normal phase in the vortex case. Lowering /i below 
some critical value \i c generates holes in the system. For 
weak disorder they form extended states, giving rise to 
a superfluid at T = (i.e. ferromagnetism for spins, 
or superfluidity for vortex systems), whereas for larger 
disorder presumably a glassy phase emerges. We will test 
these expectations by entering into the different ordered 
phases from the high chemical potential direction. In 
order to map out the phase diagram, we first consider 
the behaviour of a single hole. We then develop a many- 
body calculation for finite hole densities. 

The Mott phase is correctly described at low energies 
by the single particle characteristics. The critical proper- 
ties will be determined from the density of states (DOS), 
extending our previous formalism [[To|| . The DOS con- 
sists of a band p(Q) = —Im a(uj)/ir complete with the 
self-consistency criterion: 



a(Q) = J{{uj + hi - Jcr(w)) l )h 



(2) 



and three 6 peaks at ljsf and at i&wG with l/N 
weights. These states support superfluid and weak spin 
glass order. Their location is determined respectively by: 



where Jij = J.* a complex random number for all i < j, 
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where (■■■)h denotes averaging over the site disorder hi. 
Here u is the energy measured from the center of the 
band. Similar to the Bose glass problem [fl3 , the asymp- 
totics of the distribution of the site disorder plays a cru- 
cial role, because the first holes, being bosons, fill up the 
sites with the deepest energy levels. To investigate differ- 
ent possibilities, we generate P(h) with the asymptotics 
~ (A — \h\) r . As the DOS is the same for real J^'s, the 
subsequent results hold for XY spin glasses as well. 

For r — one finds two qualitatively different types 
of the DOS. For small values of both types of disor- 
der u>sf is distinctly below the band, supporting a su- 
perfhiid (SF) order. For strong enough bond disorder 
values the SF peak gets absorbed in the continuum, 
and the lowest energy state becomes the "Weak Glass" 
(WG) peak at uj = ujwg- ln both regimes this spin 
glass peak remains located at the lower band edge. For 
r = 1 and J = we also recover the Bose glass phase 
|jTof . For r = 2 a fourth behaviour is observed. For 
A/ J > [27/41n(4/3)] 1 / 2 = 1.393 the Weak Glass state 
is also absorbed into the band (the Supcrfluid state dis- 
solves at A/ J = 9/41n(27/16) = 1.177). This behaviour 
is present for every r > 1. In this case the self-consistent 
equation for a (Eq.|[) is satisfied by the singular terms in 
the integral over J dh P(h), which are generated by the 
deepest energy sites. In this "Strong Glass" region the 
exponent £, characterizing the asymptotics of the DOS 
p{Cu) ~ (Q — LUc)'' changes from the well-known value p"l| ] 
of Cwg = 1/2 to (sg — r - From the knowledge of the fre- 
quency dependence of the DOS one can reconstruct the 
asymptotics of the on-site Green's function at critical- 
ity. For the three glassy phases g(r) ~ t^^ +1 \ whereas 
for the superfluid g(r) decays exponentially, as the ujsf 
eigenvalue is separated by a gap from the continuum. 
The different values of these exponents again underscore 
that the transitions from the Mott Insulator into these 
four phases belong to different universality classes. 
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FIG. 1. The r — 2 phase diagram at T=0; parameters in 
units of Jo. The heavy line on the A axis is the Bose Glass. 

As the superfluid state remains separated from the 



glassy continuum by a gap, just as in the clean case, 
it is obvious that the transition into it stays mean field 
like as well jl(J . Therefore in the rest of the paper we set 
Jo = (also allowing to set J = 1) and concentrate on 
the different spin-glass transitions. We explore the differ- 
ence between these two glassy regimes by calculating the 
participation ratio P |l0| of the Weak Glass state with 
the help of the cavity method |l^] . In this approach a sin- 
gle site i is isolated, then the magnetization rrij at some 
other site j is given as the magnetization of the site in the 
absence of site i, to*, plus the effective field from i times 
the local susceptibility of rrij — to* + Uj Jji rrij, . 
Using this form in the eigenvalue equation for H, the 
magnetization is obtained self-consistently. We obtain 
for P- 1 =£X : (N/P) = {{u> WG -h-cj{Cb W G))-%. 
The participation ratio for r = 2, as plotted in the in- 
sert of Fig. 2. clearly exhibits an extended Weak Glass 
state for weak disorder, which gets localized at the same 
A c at which the Weak Glass state is absorbed into the 
continuum. 
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FIG. 2. Scaling plot of the participation ratio. 

The localization of the edge state of the continuum 
has a dramatic effect on the critical behaviour of the 
non-linear susceptibility: x^ 1 ~ 1 — (o~f)h H- As seen 
from Eq.([I|), for the Weak Glass state 1 — (of)h = 0, i.e. 
Xni appropriately diverges when \i approaches the edge 
of the band at ujwg ( m fact with a 1/2 exponent) p| . 
For the Strong Glass state Eq.(||) does not hold anymore, 
leading to the striking result that Xni is finite at this glass 
transition. 

We can also compute the ground state eigenvector of 
H numerically. Treating the participation ratio as a gen- 
eralized susceptibility suggests the scaling form: P = 
N 1 ~ v f((A — A C )A 1 / I/ ), where f(x) is a universal scaling 
function for a fixed r. The collapse of the curves for sev- 
eral system sizes in Fig. 2. is a reliable verification of this 
scaling assumption. In the fitting procedure we used the 
theoretical value of A c = [27/41n(4/3)] 1 / 2 = 1.393. The 
exponents were then determined as: r\ = 0.37 ± 0.01 and 
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v = 2.10 ± 0.10. Both exponents describe the scaling of 
spatial correlations with the finite system size. As our 
analytic technique was conceived for N = oo, presently 
the determination of their value is beyond the scope of 
this work. The one combination which can be contrasted 
with our calculation is their product, appearing in the 
form: P/N ~ (A c — A) ?,!/ . The numerical procedure 
gives rjv = 0.79, whereas analytically we get an ~ x/ In x 
asymptotics, which are numerically indistinguishable. 

So far we considered the limit of vanishing concentra- 
tion of holes in the Mott phase. Upon decreasing the 
chemical potential /i a finite concentration of holes ap- 
pears and a genuinely many-body treatment becomes 
necessary. The usual replica trick is employed, leading 
to the free energy F = lim n ^ f n /n, where: 



/„ = iTrQ 2 -l£ In *, 

i 

$i = Tr e- 0H °T T exp[a f Qa] 



(5) 
(6) 



Here the elements of the vector of annihilation opera- 
tors a are indexed (t), and that of the matrix Q 
as Qcip(t,t')- The matrix product implies a summation 
over the replica indices and an integration over the time 
variables. Also, Hf — —fa J2 a o) a a a and T T denotes the 
time ordering operator. In these equations Q satisfies the 
saddle-point condition: 



(7) 



where {•••)$. indicates averaging with respect to 

Here we introduce the key technical step of our work: 
a coherent path integral representation for the boson op- 
erators as follows: 



a 1 Qa _ 



Dip* Dip 
DetQ 



exp[— (p*Q ip + p>&* + ip*a] (8) 



With the help of these new field operators the self- 
consistent equation for Q transforms into: 



Ei(wj„) 



N ^ 



[Ei( 



(9) 



As the Fourier transforms of all matrices are diagonal 
in frequency, from now on these remain matrices in the 
replica space only. Here [g?] -1 (iw„) = (ioj n + lii)5 a p 
is the on-site hole Green's function at site i. Hi(iaj n ) 
is the self energy of the Green's function: Gi(iuj n ) a p — 
{tp* a {iuj n )tpj3(iuj n )) s it with the local action: 



Si = -y*Q V 



(10) 



We expand the action Si up to fourth order, gener- 
ating an effective tp 4 theory. Using this form one can 



construct a perturbative expansion for the Green's func- 
tion in order to determine the self energy Hi(iui n ). For 
now we stop at the level of the fully dressed Hartree dia- 
gram, which yields a replica diagonal self-energy. In the 
absence of off-diagonal elements in the self energy it is 
reasonable to proceed with a replica symmetric form for 
Q : Qapii^n) = 5 a j3a{ioj n ) + 5(iuj n )f3q. Substituting this 
form for Q in Eq.(0) determines the parameters as: 
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N ^ Ei(iuj n ) ~ a(iuj n ) 

q_ * 

N Zf{Ei{Q)-a{a)Y 



(ii) 

(12) 



We summarize the solutions of these equations at 
T = 0. In the disordered phase, characterized by q = 0, 
the self energy Ei(iu; n ) vanishes exponentially with T be- 
cause of the presence of the gap. In the dilute limit, after 
appropriate transformations, these equations are identi- 
cal with Eq.(^). 

The Weak Glass phase is easiest understood by notic- 
ing that all formulae can in fact be expanded in the on- 
site disorder strength around A = 0. In this case one 
finds E = 2, cr(0) = 1 and E(Q) ~ q ~ (ji c - fx). 
Concerning the spectrum, at low frequencies E(iui n ) = 
2 + const, x {iojn), reproducing the frequency depen- 
dence of Eq.(^|); i.e. the spectrum again starts with a 
square root form at low energies, just as at the critical 
point, as determined in the one-particle picture. Thus the 
Mott-to-Weak Glass transition is completely mean-field 
in character. In other words, the global order parameter 
q in Eq.(|l^) is generated by comparable contributions 
from all sites, indicating that the transition is driven by 
long range spatial correlations. 

The physics of the Strong Glass transition is pro- 
foundly different. As the chemical potential is lowered, 
on the right hand side (rhs.) of Eq.(|l2}), a few sites 
develop singular contributions. These local instabilities 
have to be regulated by generating a non-zero value for 
q. Obviously it is impossible to expand the rhs. around 
these singular points. This non-analyticity is well demon- 
strated in the rather non-traditional scenario for the tran- 
sition: even at the transition the coefficient of q on the 
rhs. of Eq. (|l2|) is less than one (since here Ej = /Xj, hence 
it becomes equal to {o~f)h, which is less than one out- 
side the Weak Glass phase). On the other hand, inside 
the glass phase (where q > 0) the coefficient of q must 
equal unity, therefore it jumps at the transition. This 
profound non-analyticity signals a genuinely new type of 
transition, which is driven by local instabilities. 

Next we extract the critical behaviour of q. There are 
three types of terms to the denominator Ei(0) — er(0) 
in Eq.(|^). Sites with energies larger than <r(0) give non- 
singular contributions [Xi + const, x q just as for the Weak 
Glass. Sites with energies sufficiently below cr(0) con- 
tribute ~ g 1 / 2 , whereas sites close to <r(0) yield ~ q l / z . 
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Taking into account the width of these last two integra- 
tion regions gives q ~ (// c — /i) 4 , with the subleading 
contribution ~ (/i c — /i) 6 . The spectrum is computed by 
expanding the denominator for small frequencies. One 
gets Imcr ~ cj 1 / 2 . The spectrum is quadratic in the Mott 
phase and even at the transition, i.e. the discontinuity in 
q's coefficient in Eq. (|l^) translates into a jump in the 
spectrum exponent, but the behaviour of q remains con- 
tinuous, see Fig. 3. Again, the exponents of the Strong 
Glass criticality depend on r. 
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FIG. 3. The DOS in the Strong Glass region (A = 2.0, for 
which /i c = 2.589). 

The following "phase diagram" emerges from these 
considerations. The transitions into the Weak and Strong 
Glass regions from the Mott phase are characterized by 
different exponents. The same is true approaching these 
boundaries from inside the glassy phases, however these 
two glasses do not possess different order parameters, 
hence they are not distinct phases in the thermodynamic 
sense. The spin glass order parameter q is non-zero in 
both, and the spectrum starts with a square root form in 
both as well. 

When crossing from the Weak to the Strong Glass, the 
magnitude of q is a reliable indicator, as its value drops 
by orders of magnitude (see Fig. 3). This is so because its 
asymptotic exponent of four is strongly different from the 
mean held value of one. This gives rise to an anomalously 
large critical region - the Strong Glass -, which is domi- 
nated by the local fluctuations, severely suppressing the 
value of q. When approaching from the Mott phase, the 
non-linear susceptibility does not diverge, neither does q 
develop an appreciable magnitude. Here the best indica- 
tor is the low frequency spectrum. The Mott gap shrinks 
to zero at the transition. Entering the Strong Glass the 
decreasing chemical potential shifts the spectral weight 
further down. To accomodate this weight at low posi- 
tive frequencies, the spectrum has to bulge up, taking a 
square root form: ~ (/i c — fi) x to 1 / 2 , giving way to a 
shifted quadratic form only at higher frequencies (Fig. 3). 
Thus measuring the DOS at a low fixed frequency ex- 



hibits a gap in the Mott phase, and a roughly linearly 
increasing value inside the Strong Glass. 

Lastly we comment on the issue of replica symmetry 
breaking. By explicit construction we can show that 
the replica symmetric solution is self-consistent in the 
Hartree approximation. Inspection of the structure of the 
higher order diagrams reveals that z/the replica symme- 
try is broken, the degree of this is at most proportional 
to the temperature, thus vanishing at T = 0. 

In sum, we developed a theory for the quantum vor- 
tex glass, closely related to XY magnets and bosons in 
disodered media. We included disorder both for the cou- 
pling strengths and for the site energies. Utilizing one 
particle, many body and numerical methods we deter- 
mined the phase diagram and the critical behaviour. We 
found a superfluid phase and Weak and Strong Glass re- 
gions, the physics of which is dominated by long range 
and local fluctuations, respectively. It was possible to go 
beyond previous approaches in the strongly disordered 
regime by introducing a coherent field representation for 
the hard core bosons. The most accessible experimental 
manifestation of the Strong Glass transition is the lack of 
divergence in the corresponding non-linear susceptibility. 
Analogous experimental findings have been reported in 
quantum spin-glass and vortex systems. 

We acknowledge useful discussions with S. Sachdev. 
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